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ABSTRACT 

The interstellar medium in the region of a family of cometary globules including the 
reflection nebula GM 2-39 has been analyzed basing on our observations with the slit 
spectrograph, the results of our previous observations with a Fabry-Perot interfer- 
ometer in the Ha line, Spitzer archival data together with CO microwave data. The 
structure of globules' IR emission, velocity field of ionized gas in the Ha line and of 
molecular gas in the CO line have been considered. We have detected a CO cavern 
around the eastern globules and faint high- velocity Ha features of surrounding gas. 
The most probable sources of ionizing radiation and wind of the Cyg OBI associa- 
tion responsible for the globules' formation are proposed. Based on our multicolour 
photometry, it has been found that the brightest compact source in the southern 
'head' of the cometary globule - the star HBHA 3703-01 (IRAS 20153+3850) - is a 
B(5-6)V star with E(B — V)=1.18 mag. The strong and broad Ha emission line was 
revealed in its spectrum. The spectral energy distribution of HBHA 3703-01 in the 
0.44-24 /im range has been modeled. It is shown that the star has a hot dust envelope 
with Tdust = 1400 K and ry = 1.1. Besides, the star HBHA 3703-01 illuminates the 
encompassing diffuse nebula GM 2-39 with a diameter of about 30 arcsec. With regard 
to all obtained observational data for HBHA 3703-01, this object may be classified as 
the Herbig Ae/Be star. 

Key words: ISM: kinematics and dynamics - ISM: clouds - ISM: lines and bands - 
infrared: ISM - ISM: individual objects: GM 2-39 - stars: individual: HBHA 3703-01 



1 INTRODUCTION 

The concept of a globule was introduced as early as in 1940s 
for 'dark nebulae' - opaque compact clouds with sharply de- 
lineated boundaries on the backgrou nd of emission n e bulae . 
Desc ribed for t he fir st time by lHawarden &; Brand! (| 19761 ) 
and lSandauistl (|l976h , cometary globules (CGs) are charac- 
terised by cometlike morphology - a compact opaque dust 
'head', frequently surrounded by a bright rim, and a long 
dust 'tail'. It is supposed that such a shape of CGs is formed 
as a result of the action of ultraviolet radiation flReipurthl 
Il983h and /or of shock waves initiated by supernovae or stel- 
lar winds (|Brandlll98ll : lBrand et aHll983l : ISaim et al.lll988l : 
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Pitt ardll201lh on small dense clouds immersed in a rarefied 
medium. As a rule, CGs' tails are oriented radially from the 
source that has formed them. The problem of the origin of 
CGs is still discussed. The study of their nature is of special 
interest because of their connection with the process of star 
formation: in heads of many globules compact IR sources - 
young stellar objects (YSOs) - have been detected. 

The launch of IR space observatories opened a new 
stage in the study of CGs formation processes and their 
generic connection with star-forming regions. Dense dust 
clouds and their 'fine structure', formed by radiation and 
wind from young stars and determining the process of for- 
mation of a new generation of stars, are studied today very 
intensely using IR observations. The most impressive results 
have been obtained for the Carina Nebula (NGC 3372) from 
observations with t he Spitzer and Hubble space telescopes 
([Smith et al.ll2010al lbh. Here, multiple new Herbig-Haro ob- 
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jects with distinct bipolar jets, cometary globules with signs 
of star formation, circumstellar gas-dust discs (proplyds) 
have been detected. The spatial distribution of IR emission 
in different spectral bands revealed a clear picture of con- 
tinuing star formation in dust pillars, which result from the 
action of young stars on the remnants of the parent molec- 
ular cloud. 

In this paper we study the structure and kinematics of 
the interstellar medium towards a family of CGs (CK2000 ~ 
20 h 15 m 40 s -20 h 17 m 45 s , feooo ~ 38°57 / -39°10 / ) detected in 
the archival data of the Spitzer Space Observatory in the 
field of the supershell in Cygnus. We carried out spectral 
observations of a compact stellar object in the head of the 
southern globule. The purpose of this work is the analysis of 
kinematics of ionized and molecular gas, revealing a physical 
association between different components of the stellar and 
gas populations of the region and search for traces of the 
Cyg OBI stars' wind action on cometary globules. 

The description of IR images of CGs and their compar- 
ison with optical images is given in Section [2] In Section [3] 
we report the results of our study of the kinematics of the 
ionized gas towards the family of CGs using the data of our 
long-term interferometric Ha observations of the gas-dust 
complex in Cygnus. In Section [4] we present evidence for the 
association of the family of IR globules with a complex of 
molecular clouds and stars at the periphery of Cyg OBI; 
we have searched for probable kinematic signs testifying to 
the action of stellar wind on gas near CGs. In Section [5] we 
present the results of photometric and spectral observations 
of the star HBHA 3703-01 and the model of the dust enve- 
lope around this star. In Section [6] we summarize the main 
results together with their implications. 



resolution and, probably, to the fact that heated dust re- 
sponsible for the main emission in this band can be present 
also in the neighbourhood of the globules.) The IR images 
of the globules coincide with the regions of opaque dust in 
Fig. [TJl surrounded with gas emitting in the optical range. 
The bright contours of the IR globules are inside their op- 
tical boundaries. 



In the IR range as in the optical the 'three-head' struc- 
ture and the tail of the eastern IR globule are visible; the 
size of the largest head reaches 3 arcmin. The heads are 
'filled' with IR emission of different structure: the north- 
ern oval head is limited with a bright contour and has 
a cellular IR structure (Fig. [TJ)); in the 'middle', cone- 
shaped one, the brightness of IR emission decreases to- 
ward the tail. The brightest southern head hosts the re- 
flection nebu la GM 2-3 9 detected in the optical range by 
Gyulbudaghian & Maghakian (1977) (Fig.[T]l). Its size does 
not exceed 1 arcmin, and the integrated stellar magnitude 
on Palomar Sky Survey images is nearly identical in the red 
and blue light. GM 2-39 is the brightest extended IR source 
in the area; its morphology in the optical and IR ranges 
is identical (Fig. [TJ)). In the western IR globules having an 
elongated shape and smaller size, we reveal bright, some- 
times intermittent contours with a brightening in the region 
of the supposed head and at the end of the tail. (Note that 
the studied objects are not classical cometary globules, typ- 
ical of which is not only a distinct bright head, but also a 
long dust tail. Probably, the western globules are fragments 
of the tails of the eastern heads; a 'classical' cometary nebula 
is only the southern part of the eastern globule.) 



2 FAMILY OF COMETARY GLOBULES IN 
THE INFRARED AND OPTICAL RANGES 

The family of cometary globules in the region / ~ 
76°30 / -76°40 / ; b ~ 1°55 / -2°15 / was detected in 
the infrared using observational data of the Spitzer 
Space Observatory obtained from th e Spitzer L egacy 
Survey of the Cygnus-X Comple£| (JHora et al.1 I2008L 
http : //www . cf a . harvard . edu/cygnusX) . The archival data 
were downloaded from the Spitzer Heritage Archival. We 
had at our disposal images of the region in five wavebands: 
3.6 jiim, 4.5 /mi, 5.8 /mi, 8.0 /mi observed with the InfraRed 
Array Camera (IRAC) and 24 /xm with the Multi-Band 
Imaging Photometer for Spitzer (MIPS). The spatial reso- 
lution of the images obtained with these cameras is 0.6 and 
2.45 arcsec per pixel, respectively. We composed mosaics of 
the images using the MOPEX softwarqj. 

Figures [TJibc show the images of the region in the 3.6, 
8 and 24 /mi bands. It is visible that the family of IR glob- 
ules consists of several lined in a chain, identically oriented 
CGs with sizes from ^ 1 to 7 arcmin. The morphology of 
the IR globules in all the wavebands is the same as a whole. 
(The 'blurred' image in the 24 /xm band is due to a lower 
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Peculiarities of the globules in the IR range are thin 
small filaments in the form of 'whiskers' near the heads 
(Fig. [TJib). Similar IR filaments have been detected in a 
large number of cometary globules. As an exam ple, let u s 
mention globules in the Carina Nebula (Smi th et al.ll2002n . 
The filaments in front of the southern head are curled to- 
ward the tail; they resemble rims noticeable in the optical 
range at the peripherie s of the he ads of cometary globules in 
the Gum Nebula (|Reipurthlll983h . From the east the 'Three- 
headed' IR globule is also bordered with an optical rim, a 
much broader one than that of IR (Fig. [TJi). Between the IR 
head and its optical rim a sharp boundary, especially promi- 
nent in the case of the southern head, is seen (Fig.[T|3). The 
western globules are limited with bright IR contours; how- 
ever, there are no optical counterparts for them. Note that 
there are regions of cool dust in the tails of the eastern and 
one of the western globules (Fig. [1]). 



Thus, the morphology of CGs in the IR bands is as a 
whole consistent with the optical images: dust head and tail, 
restricted with bright IR and/or optical rims. The optical 
emission encircles the IR globules. However, in the IR range 
there is a characteristic feature that has no optical coun- 
terpart: thin filaments in front of the heads of the 'Three- 
headed' CG. The structure of other globules in the IR range 
is also more expressive than in the optical, since at the pe- 
riphery the globules are bordered with bright IR rims in all 
considered spectral bands. 



A family of cometary globules 3 



39° 08' 








Pfe^^V 


39° 00' 




L ^HI^SBBSlSsv 







39° 08' 



39° 00' 



39° 08' 



39° 00' 



20 h 17 m 00 



20 h 16 m 00 




20 h 17 m 00 s 



20 h 16 m 00 s 




39° 08' 




20 h 17 m 00 






20 h 16 m 00 s 









k 'I." "»'•". , •■' ' ■-" • '• • • »„v,.'/ .'■ ' . . . '-"aP V * • ■ ' " 

■■■■'• i-5 ♦ :,;-/■: /- :■•; • ■'•' -V '. ■•"■ . • -. ■--■ ■—.■::.;:/". . * •;,;;■ 

V.*# "-"V-: ■■■ >- '■:,..■■ ;.■.>;.■ v> ■■".■■ o : .£ " 

"•• .'• ■■**»• *Av. .' . ■',"•••'•* •'.":' * * \» ." . ' '-' ** • -. . 

SSto'! li^^ ;WI«i 



20 h 17 m 00 s 



20 h 16 m 00 s 



Figure 1. The family of cometary globules in the bands 3.6 /im (a), 8 /im (b), 24 /im (c) according to the Spitzer data (the IR image 
shown by red is combined with the optical shown by blue) and in the optical range (Palomar Sky Atlas, E-map) (d). Right lower corner 
of Fig. b: enlarged map of the southern head of the eastern globule with the embedded cometary nebula GM 2-39; a cross marks the star 
HBHA 3703-01. The localization of the slit spectrogram is indicated (Fig. d). 
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Table 1. Radial velocities of ionized hydrogen 



Location 



(i) 



l,b Nl ((V LS r) - A) . . . N2 

02000,(52000 ((Vlsr) + A), km s _1 

(2) (3) (4) (5) 



Eastern globules 



Western globules 



76°32 / -76°40 / 36 

l°52 / -2°06 / 



20 h 15 m 58 s -20 h 16 m 50 s 
39°05 / -39°09 / 



30 



CO cavern around WR 139 20 h 17 m 00 s -20 h 20 m 00 s 230 
(region II) SS^O'-SS'MS' 

CO cavern near HD 193595 20 h 18 m 18 s -20 h 19 m 20 s 63 
and HD 228841 (region III) 38°52 / -39°10 / 
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3 VELOCITIES OF IONIZED GAS IN THE 
REGION STUDIED 

Detecting physical connection between different components 
of ionized gas, dust, dense molecular clouds and stellar pop- 
ulation of the extended region around the considered CGs 
requires their radial velocities. It is difficult to find such a 
connection, because the region belongs to a heavily popu- 
lated field, that is viewed along the spiral arm. 

In the analysis of velocities of ionized gas towards CGs 
and their neighbourhood we used the results of our observa- 
tions of the Cygnus gas-dust complex in the Ha line carried 
out with a Fabry-Perot interferometer (IFP) mounted in the 
Cassegrain focus of the 125-cm reflector of the Crimean Lab- 
oratory (Sternberg Astronomical Institute, Moscow State 
University) . A technique of o bservations and their pr ocess- 
ing were described in detail bv lLozinskava et al.l (|l998h . The 
field of view and angular resolution of the observations was 
10 arcmin and 3-4 arcsec, respectively; the actual spectral 
resolution was 10-15 km s _1 . The line profile was approx- 
imated with one or several gaussian curves in the supposi- 
tion that the FWHM of each component is larger than the 
halfwidth of the instrumental profile and signal-to-noise ra- 
tio is ^ 5. 

In 1993-1995 we obtained 15 IFP-images in a broad 
vicinity of CGs. Here, as everywhere in the Cygnus com- 
plex, multipeak profiles of the Ha line are observed. They 
consist of one or several bright main components with ve- 
locities of the peaks within the interval Vlsr ~ 0-20 km s _1 
and of faint 'shifted' components in the line wings with ve- 
locities up to Vlsr ~ —80 km s _1 and Vlsr ~ 60 km s _1 . 
The radial velocities of the main component are due to sys- 
tematic motions created by the galactic rotation and spiral 
density shock waves. Those of the shifted components are 
caused by local high- velocity flows of gas related to shock 
waves initiated by supernovae and/or stellar winds. At large 
negative velocities we cannot eliminate also a contribution 
from distant HII regions of the Perseus Arm. 

For allocation of characteristic systematic motions of 
ionized hydrogen in the field of CG and their neighbour- 
hoods we have plotted 'diagrams of occurrence' of radial 
velocities in selected areas. Using them, we have found ra- 



dial velocities of gas Vlsr characteristic of the region as 
the average values at half-maximum of the corresponding 
diagram together with root-mean-square errors of their es- 
timate A. The obtained velocities are listed in Table [1] in 
columns: (1) region identification (CO caverns are shown in 
Fig. [2j)); (2) Galactic or equatorial coordinates; (3) num- 
ber of points Nl at which radial velocities were measured; 

(4) range of the prevalent radial velocities of systematic mo- 
tions ((Vlsr)— A) • • • ((Vlsr) + A) (emphasized in boldface); 

(5) number of measurements N2 from which the velocity 
range was determined. Given for each of the four areas in 
the first line are velocities of systematic motions of gas, in 
the second and third lines are velocities of shifted positive 
and negative components of the Ka line. The field of radial 
velocities in the broad neig hbourhood of globule s con sidered 
here w as analyzed also by lEsipov et al.l (| 19961 ) and ISitnikl 
(|201lh . 

The detected intervals of systematic radial velocities 
of ionized hydrogen are typical of the population of the 
Cygnus Arm and are consistent with the corresponding ve- 
lo cities of molecular clouds of this arm according to the data 
of iLeung &; Thaddeusl ()1992| ) (see also ISitnik &; Lo zinskava 
2009) and with modern CO observa tions of the Cyg X re- 
gion with a high angular resolution (jGottschalk et al.ll2012l 
and references therein). 



STELLAR POPULATION, STRUCTURE 
AND KINEMATICS OF GAS IN THE 
VICINITY OF THE COMETARY GLOBULES 



According to adopted concept s (Reipurth 1983^ Brandlll98ll : 



iBrand et aHll98alSahu et al.|[l988l : IPittardll201lh . cometarv 
globules are shaped by shock waves and/or ultraviolet radi- 
ation, so the head should point to the energy source, and 
the tail is extended radially from it. In our case individual 
CGs as well as the entire family of CGs lined in a chain are 
oriented such that the heads 'look' towards the plane of the 
Galaxy (see Fig. [2ji). Accordingly, the energy source should 
be sought for in a region southeast/east from CGs. Note 
that the considered region includes many optical globules 
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Figure 2. Neighbourhood of the family of cometary globules, (a) - The optical map in the red range. Numbers 1, 2 and 3 label the 
stars HD 228841, HD 193595 and WR 139. Rectangles limit the areas for which the P/V diagrams were plotted. The nebula Simeiz 55 
is in the area (0=2000 ~ 20 h 13 m - 20 h 18 m , £ 2 qoo ~ 38°30 / - 39 o 20M (b) - Image of the region in the 24 /im band with contours of the 
CO emission at Vlsr — 2.5 — 5 km s -1 , data of iLeung fe Thaddeusl Jl992). Contours are labeled in the order of increasing brightness: 1 
corresponds to the weakest CO emission. Roman numerals: CO caverns mentioned in the text and Table [T] Ticks mark the same stars 
as in Fig. a. 



(Fig. [T]i, [2^l); however, not all of them show up in the IR 
range and not all are surrounded by bright optical rims. 

We can reveal toward the globules and their neighbour- 
hood in the southeast sector of the complex the following 
components of the interstellar medium and stellar popula- 
tion. 

1. A CO cavern (Z - 76°30'- 76°40', b ~ 1°55'- 2°15 / ) 
i s observed at the velocity Vlsr — 2.5—7.5 km s _1 
([Leung &; Thaddeusl Il992h towards the eastern 
globules (region I in Fig. [2)3). According to Table [1] ionized 
hydrogen that the globules are embedded in radiates at 
the same systematic velocities Vlsr — 4-6 km s _1 . In 

a 36 arcmin region including the family of globules, the 
emission in the HI 66a radio recomb ination line at v elocity 
Vlsr — 5.4 km s _1 is also observed (|Landeckerlll983h . 

2. CGs are at the periphery of the young 
stellar association Cyg OBI, inside a large-scale 
shell around Cyg OBI and Cyg OB 3 detec t ed in 
the optical, radio and IR ranges (Brand fe Zealev 19771 ; 
iLozinskava fe Sitnild Il988l : ILozinskava fe RepiiJll99Ch. The 
distance to Cyg OBI i s 1.5 kpc (|Garmanv fe Stencellll992l ; 
ISitnik fe Mer"niklll99d ). The median radial velocity of the 
stars of Cyg OBI found from the velocities for 34 of 70 stars 
of the associa tion is 4 km s _1 with a dispersion of 8.9 km s _1 
(Sitnik 2003), which is consistent with the quoted above ve- 
locities of ionized gas and CO emission and allows us to 
suggest a close spatial location of Cyg OBI and globules. 
Note that in the studied direction of the Cygnus Arm, ex- 
cept for the cluster Be 86 which is a member of Cyg OBI 
([Forbes et alJll992T h no other young stellar group has been 
detected. 

3. At a separation of 15 — 30 arcmin to the east 
and southeast from the chain of CGs, on its ex- 
tension, at least three stars with a powerful wind 



are localized: HD 228841 (Of), HD 193595 (Of) and 
WR 139 ( Fig. Hi). Both Of stars are members of Cyg OBI 
association (Blaha fe Humphrevd l 1989k iGarmanv fe Stence 
1992), and WR 139 is its probable member ([Van der Huch 
l200ll ). All three sta rs are observed in a reg ion with de- 
pleted CO emission ([Leung fe Thaddeusl 1 1992h detected at 
velocities Vlsr > 2.5 km s _1 inside the supershell around 
Cyg OBI, Cyg OB3 (regions II, III in Fig. [2]o). This 
CO cavern, which is probably formed by combined action 
of stellar wind and radia ti on of the s tars of Cyg OBI 
(|Lozinskava fe Sitnild fl988l : ISitnikl [JOllh . is limited from 
the north by a complex of molecular clouds at velocity 
Vlsr ^ 2.5-7.5 km s _1 . The star HD 228841 is on the exten- 
sion of the chain of CGs and is separated from the nearest 
heads of CGs by 6 pc (see Fig. [2]). The star WR 139 is lo- 
cated in the same direction at a projected distance of 13 pc 
from the heads of CGs. The star HD 193595 is located at 
9 pc to the east from the family of globules; the chain of 
CGs is partially deployed toward it. The extension of the 
chain of CGs at the distance of Cyg OBI is 9 pc, and the 
sizes of the globules are not larger than 3 pc. 

4. Toward the family of globules high- velocity mo- 
tions of ionized hydrogen with characteristic velocities 

Vlsr 56... - 30 km s" 1 and Vlsr ~ 33-44 km s" 1 

are observed (see Table [J). Earlier, high- velocity motions of 
ionized hydrogen in this re gion were attribute d to the fila- 
mentary nebula Simeiz 55 (|Esipov et alJll99q ) (Fig. [2^). In 
the sky plane the nebula is a part of the supershell around 
the Cyg OBI association. The family of globules is observed 
on the background of this nebula (Fig. [2^l). The nature of 
high- velocity motions and filamentary morphology of Simeiz 
55 are not completely clear, because in the nebula neither 
synchrotron radio emission, which would testify to a super- 
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Figure 3. P/V diagrams to the northwest from HD 228841 (a) 
and to the west from HD 193595 (b) for the regions designated 
in Fig. [2jt- Dashed lines denote borders of the regions I and III 
shown in Fig. [2j). 



nova explosion, nor a star with a powerful wind have been 
detected. 

High- velocity motions of ionized hydrogen in the 

neighbourhood of the family of globules could be evidence 
for the physical association of cometary globules with nearby 
stars of the Cyg OBI association. To find probable kinematic 
signs of the Of stars wind action on gas in the neighbour- 
hood of the globules, we have analyzed high- velocity com- 
ponents in the Ha line. Along the chain of the globules, we 
have plotted distributions of radial velocities as functions of 
projected distances to each of the stars (so-called P/V dia- 
grams) (Fig. [3]). The regions for which P/V diagrams were 
plotted are shown in Fig. [2^i. 

In Fig. [3^l high velocities of the shifted components in 
the Ha line (up to —60 .. .—80 km s _1 ) are found at a separa- 
tion of 0?1 (2.5 pc) from the star HD 228841 and are visible 
up to the supershell around Cyg OBI, where they become 
equal to the mean velocity of the surrounding medium. We 
see no definite decrease in the radial velocity that could be 
attributed to a head-on collision of the shock driven by the 
wind of HD 228841 with the heads of the globules. Indirect 
evidence for the interaction between the wind of the star 
and the family of CGs can be the shape of the heads pro- 
late perpendicularly to the direction toward the star as well 
as the presence of the CO cavern surrounding the eastern 



globules (see item 1) and of fine IR filaments curled away 
from the star. The existence of high- velocity motions ob- 
served beyond 0?25 (7 pc), that is farther than the heads, 
can be explained by gas circumflowing the family of CGs. 
The main non-shifted component of the line in Fig. [3] sums 
up the emission of non-perturbed gas near the star and in 
the line of sight. 

In Fig. [3j) high velocities, up to 60-80 km s _1 , are 
visible at a projected separation of 0?45 (12 pc) from 
HD 193595, and decrease by the absolute value to the veloc- 
ity of the surrounding medium at 0?8 (21 pc). The family of 
globules is within the limits of this distance. Note also that 
HD 193595 and the family of CGs are at the opposite sides 
of the CO cavern detected at velocity Vlsr > 2.5 km s _1 
(region III in Fig. [2j)). In Fig. [3}) this rarefied region of the 
interstellar medium extends to not less than 0?25 (7 pc) 
from HD 193595. Therefore, high velocities of ionized hy- 
drogen are found only far from the star. 

Thus, the family of CGs falls within the zone of action 
of the wind from the Of stars HD 228841 and HD 193595. 
The surrounding ionized gas swept by the wind of these 
stars expands at a velocity of up to 60-80 km s _1 and cir- 
cumflows the family of cometary globules, since we see high 
velocities along the entire chain of CGs. We have found no 
obvious kinematic signs of braking of high- velocity gas due 
to interaction with the Three-headed globule. 

From what was said it follows that the generic associ- 
ation between the family of CGs and the above-mentioned 
stars of Cyg OBI is possible. This is based on the following 
arguments. The family of CGs is surrounded by ionized gas 
radiating at systematic velocities Vlsr — 3-7 km s _1 . The 
median radial velocity of the Cyg OBI stars is 4 km s _1 . 
At the same velocities a CO cavern is observed around the 
eastern globules (region I in Fig.[2]3). This small cavern was 
formed in a complex of clouds limiting another CO cavern 
detected at Vlsr > 2.5 km s _1 and connected with the 
Cyg OBI association (see regions II and III in Fig. [2}) and 
in Table [T}. The above-mentioned Of and WR stars belong 
to Cyg OBI. In the sky plane the stars lie on the exten- 
sion of the chain of CGs within the limits of the action of 
their UV emission and wind. Important signs of the prob- 
able physical association of these stars with the cometary 
globules are radial elongation of the chain of CGs relative to 
the stars and its morphology: the heads with a bright rim 
point toward the stars, large bright globules are closer to 
the stars than small ones. IR rims observed near the heads 
of globules and cu rled toward the ta ils belong to type 1 in 
the classification of ISmith et al.l (|2002[ ) (see Fig. 33 therein). 
This is one more argument for the association of the stars 
and globules, since such a shape is formed under the action 
of wind and radiation of a massive star or cluster. 



5 THE STAR HBHA 3703-01 AND THE 
NEBULA GM 2-39 ASSOCIATED WITH 
COMETARY GLOBULES 

The nebula GM 2-39 hosts the star GSC 03151-00990 
(a = 20 h 17 m 08?05; S = 38°59'29"4 (2000)) identified with 
the infrared source IRAS 20153+3850 (Fig. lbd). The star 
i s also included in the cata log of Ha emission objects 
(JKohoutek fe Wehmeverlll999h as HBHA 3703-01. Accord- 
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Figure 4. The image in the R band of HBHA 3703-01 with com- 
parison stars. 



ing to the classification of emission spectra of stars in this 
catalog, the Ha emission line in the spectrum of HBHA 
3703-01 has a medium intensity on the background of a 
moderate continuum. The catalog lists the star' brightness 
of V = 13.6 mag. 

To get more complete information about the star and 
surrounding interstellar medium we carried out the spectral 
and photometric observations. 

5.1 Photometric observations of HBHA 3703-01 

The photometric observations were carried out with the 60- 
cm reflector of the Crimean Observatory of the Sternberg 
Astronomical Institute equipped with Apogee AP-47p cam- 
era. A total of 158 images in the UBVR C RI bands were 
obtained in September-November 2012. All reductions were 
made using IRAFq MAXIM and programs of V.P. Goran- 
skij. 

The comparison stars were calibrated with respect to 
the star HD 228743. Its stellar magnitudes (U = 11.71 mag, 
B = 11.51 mag, V = 11.22 mag, R c = 11.04 mag, 
R = 10.93 mag and I = 10.72 mag) were obtained us- 
ing the following pho tometric standards : PG2 336, PG0231, 
RU149 riLandoltlll992h M92, NGC 7790 (|Stetsonl2000h . M81 
([Richmond et al.l ll994). The red uction o f R c to R and I c to 
I was done with equations from IBesselll (| 19791 ). The image 
of HBHA 3703-01 and comparison stars are shown in Fig. [4] 

The colour terms for transformation of instrumental 
magnitudes ubvr c i to standard UBVR C I were determined 
using observa tions of standard s in M 92 and NGC 7790 
(jStetsonl IJOOoh and in M82 ([Richmond et alJ Il994h . The 
equations u = U + K U (U-B) + C U ] b = B + K b (B-V) + C b ; 
v = V + K V (B -V) + C v ; r c = R c + K rc (V - R) c + C rc 
and i = I + Ki(R c — I) + C% were solved for colour terms: 
K u = -0.08, K b = 0.04, K v = -0.03, K rc = -0.04 
Ki = 0.027. 

The results of our photometric observations of HBHA 
3703-01 are presented in Table [2] The star brightness varies 



4 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by AURA under cooperative agreement 
with the National Science Foundation 
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Figure 5. HBHA 3703-01 in the U - B, B-V diagram The 
solid curve is the main sequence according to Straizvs (1977); the 
dashed line represents the reddening line in form E(U — B) / E(B — 
V) = 0.796 + 0.020£(B - V). 



from night to night, but these variations do not exceed 
3<r; therefore, the conclusion on star photometric variabil- 
ity would be premature. 

The comparison of mean colour indices of HBHA 3703- 
01, U - B = 0.4 ± 0.14 mag, B - V = 1.01 ± 0.13 mag, 
with normal colour indices of main-sequence stars shows 
that HBHA 3703-01 is a B(5-6)V star with the colour excess 
E(B—V) = 1 . 18=b0. 17 mag. For that we used intrinsic colour 
indices from IStraizvsl (|1977l ) and the reddening line of Cyg 
OB2 in for m E(U - B) /E( B - V) = 0.796 + 0.020^(5 -V), 
derived by I Turner! ()1989l ). We note that there are only 
three stars HD 193032 (Bl III), HD 192422 (B0.5I) and 
HD 228841 (07.51) with known UBV-daia and the two- 
dimensional spectral classification in a close (r ^ 30 x ) vicin- 
ity of HBHA 3703-01. Photometry of these stars provides 
an average E(U - B)/E(B -V) = 0.81 ± 0.01 that is in a 
good agreement with Turner's data. The position of HBHA 
3703-01 on the U - B, B - V diagram is indicated in Fig. El 

We estimated the visual absolute magnitude of HBHA 
3703-01 using our data -V= 14.98 ± 0.07 mag and E(B - 
V) — 1.18 =b 0.17 mag - assuming that the star belongs to 
the Cyg OBI complex at the distance 1.5 kpc. It turned out 
to be My — 0.4 ± 0.1 mag, in a good agreement with the 
visual luminosity of B6 s tars on the zer o-age main sequence 
(ZAMS, M v = 0.3 mag. lStraizvslll982h . 

The colour excess, E(B — V) = 1.18 mag, derived by 
us, characterises the total light extinction involving the in- 
terstellar and circumstellar components. The problem of 
light extinction toward HBHA 3703-01 is complicated, since 
the star is embedded in an extremely inhomogeneous dust 
medium. The nearest to HBHA 3703-01 star with the known 
reddening- HD 228841 (Of7.5p) with E(B-V) = 0.87 mag, 
member of Cyg OBI - is located at about 18 arcmin from the 
emission star outside the region of large extinction. Adopt- 
ing for HBHA 3703-01 this value as its interstellar extinc- 
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Table 2. UBVR C RI magnitudes of HBHA 3703-01 from observations in September-November 2012 



JD-2400000 



U 



°U 



B 



°B 



V 



°v 



R c 



a Rc 



R 



°R 



56187.3 
56216.2 
56218.2 
56219.2 
56240.2 
56243.2 
56247.2 



07 



15.89 


0.05 


14.96 


0.01 


- 


- 


- 


- 


- 


- 


15.99 


0.04 


14.99 


0.01 


- 


- 


13.89 


0.01 


13.22 


0.01 


15.95 


0.06 


14.89 


0.02 


- 


- 


- 


- 


- 


- 


16.02 


0.03 


15.01 


0.02 


- 


- 


13.86 


0.01 


13.26 


0.02 


16.02 


0.02 


14.94 


0.01 


14.29 


0.01 


13.87 


0.01 


13.29 


0.02 


16.06 


0.02 


15.03 


0.03 


14.29 


0.01 


- 


- 


13.27 


0.02 


15.97 


0.02 


15.00 


0.02 


14.28 


0.01 


- 


- 


13.17 


0.02 



56252.1 


16.45 


0.08 


16.04 


0.03 


15.02 


0.05 


14.27 


0.01 


- 


- 


13.30 


0.02 


mean 


16.45 


0.08 


15.99 


0.10 


14.98 


0.07 


14.28 


0.02 


13.87 


0.02 


13.25 


0.04 



tion, the circumstellar reddening component of HBHA 3703- 
01 may be estimated as about 0.3 mag. 

The reflection nebula GM 2-39 surrounding the star 
HBHA 3703-01 is well visible on images in all bands. The 
measured colour indices of the nebula, B—V — 0.5±0.2 mag, 
V — i? = 0.7 ± 0.2 mag, are considerably bluer than that of 
the star: B-V = 1.01 ±0.13 mag, V-R= 1.11 ±0.11 mag, 
due to the strong scattering of the stellar light in this nebula. 

We note the star H BHA 3703-01 was too much brighter 
in lDroege etal] <l2007h (V = 13.4 4 ± 0.12 mag) than in our 
observations. iDroege et al] (|20070 made the measurements 
with a 30" aperture so their data refer to the total star's and 
a considerable part of reflection nebula GM 2-39 radiation. 



5.2 Spectral observations of HBHA 3703-01 and 
of its environment 

Long-slit spectra at low spectral resolution were obtained 
with a large- aperture spectrograph in the Cassegrain focus 
of the 125-cm reflector of the Crimean Laboratory of the 
Sternberg Astronomical Institute in 2011 and 2012 years. 
We used a diffraction spectrograph with a 600 lines/mm 
grating. The observations were conducted with a slit about 
8 arcmin long and 4 arcsec wide; the scale along the slit was 
1.98 arcsec per pixel. As a detector an ST-402 camera was 
used (array size 765x510 pixels, pixel size 9x9 /mi). The 
actual spectral resolution was FWHM = 6.0 A. The typ- 
ical achieved signal-to-noise ratio (S/N) is about 10. The 
spectrograph slit was directed east-west (the localization of 
the spectrograms is shown in Fig. Id). Table [3] gives the log 
of observations: for each spectrum it lists the date, spectral 
range AA and total exposure Te XP . The spectrophotom etric 
standards 40 Cyg and 57 Cyg ([Glushneva et al.l 11998) were 
used for flux calibration. The spectra were processed by the 
standard CCDOPS software as well as by SPE code cre- 
ated by S.G. Sergeev in the Crimea n Astrophysical Obser- 
vatory (jSergeev &; Heisbergerlll993T ). The primary process- 
ing included dark current subtraction, extracting individual 
spectra from CCD images and their wavelength calibration. 
The star HBHA 3703-01 was placed at the slit centre. 
The only line detected in the observed spectrum of the star 
is broad Ha (FWHM = 9.5 ± 0.5 A). After excluding the 
instrumental profile, FWHM of the line is 7.3 ± 0.5 A; this 
corresponds to a high velocity, Av = 169 ± 25 km s _1 . The 
large linewidth can be caused by stellar wind, star rota- 
tion or turbulent motions in the outer envelope of the star. 
We measured the equivalent widths of Ha on four spectro- 



Table 3. Log of spectral observations 



Date 



JD-2400000 



AA, A 



^exp; > 



05/06.10.11 


55840 


26/27.10.11 


55861 


15/16.09.12 


56186 


10/11.10.12 


56211 


15/16.10.12 


56216 



5700-7200 600 

5700-7200 900 

4000-5700 1800 

4000-5700, 5700-7200 1800, 1800 

4000-5700, 5700-7200 1800, 1800 



grams obtained in 2011 and 2012 and did not find any vari- 
ations exceeding the measurement errors. The mean equiva- 
lent width is W(Ra) = 103 ± 2 A. Most likely, the emission 
line is formed by the collisional excitation mechanism. By 
this, the object is similar to Herbig Ae/Be stars, character- 
istic of star- forming regions. 

We measured also the absolute flux in Ha for HBHA 
3703-01 F(Ha) = (5.1 ± 0.5) x 10" 13 erg cm'V 1 . 
When corrected for light extinction with Ay — 3.1E(B — 
y)=3.65 mag, we find F (Ha) = 7.0 x 10" 12 erg cm" 2 s _1 . 
Adopting distance to the star d = 1.5 kpc found for the 
Cyg OBI association, we obtain the luminosity in the line 
L(Ha) = 47rd 2 Fo(HaO = 1.8 x 10 33 erg s _1 , or 5.9 x 10 44 pho- 
tons per second. A typical B(5-6)V star delivers much less 
photons in the range A < 912 A; therefore, the Ha line can 
not be produced by a process of photoionization followed 
by recombination. As will be shown below, HBHA 3703-01 
may be classified as the Herbig Ae/Be star. The origin of 
the emission lines in such stars is still under discussion. 

In the spectrum of ionized gas, along the entire spectro- 
graph slit (~8 arcmin), the emission lines H/3, [Olll] A 4959, 
5007 A, Ha, [Nil] A 6548, 6584 A, [Sn] A 6717, 6731 A are 
seen. Figure [6] shows the spectrum of the star HBHA 3703-01 
after subtraction of the nebula background and the emission 
spectrum of ionized gas. 

The reflection nebula GM 2-39 has a diameter of about 
30 arcsec and is well seen in the optical and infrared, but 
it is invisible on our spectrograms owing to its low surface 
brightness. 

The emission spectrum of ionized gas of the surround- 
ings is due by its origin to the general field of ultraviolet ra- 
diation of the hot stars WR 139, HD 193595 and HD 228841. 
The intensity of the emissions vary along the slit with the 
inhomogeneity of interstellar medium. Figure shows the 
distribution of the Ha flux of ionized gas along the slit from 
east to west. At 1.3 arcmin to the east from the star HBHA 
3703-01 there is a dense optical filament = rim of the south- 
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Figure 6. Spectrum of the star HBHA 3703-01 (black line) and emission spectrum of ionized gas in the region of the filament bordering 
the southern head (red line). 
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ern head of the globule. On the slit at this place the emis- 
sion lines of hydrogen, sulfur and nitrogen are especially 
enhanced. 

We measured the intensity ratio of basic emission lines 
of ambient ionized gas. The average values along the slit are: 
F(5007)/F(H/3) = 1.15 ± 0.17, F(6584)/F(Ha) = 0.31 ± 
0.04, F(6717 + 6731)/F(Ha) = 0.23±0.02. These values are 
typical of the interstellar medium with solar abundances and 
they fall into the position of typical H n regions on diagnostic 
diagrams log(|Nll l/Ha) - log([Olll1/H ff) and log([Sll]/Ha) 
-log([Om]/H/3) (|Kniazevet al.ll2008h . 

The mean value of F(A6717)/(A6731) = 1.47 ± 0.10 is 
very close to its limiting value for low densities, allowing to 
determine only an upper limit of the electron density in the 
ambient gas, N e ^ 40 cm -3 . 



5.3 Model of the spectral energy distribution of 
HBHA 3703-01 with a dust envelope 

The object HBHA 3703-01 is well visible on original Spitzer 
images in the 3.6, 4.5, 5.6, 8, 24 and 70 /mi IR bands. 
On all images except for 70 /xm it looks as a starlike ob- 



ject surrounded by a nebula. On the last image only the 
nebula is visible. Using these images except for 70 /xm, 
we estimated its fluxes in each of the bands in units of 
10~ 16 erg cnrV 1 !- 1 : F 3 . 6 = 39.2, F4.5 = 28.8, F 5 . 6 = 
29.1, F 8 = 17.1 and F 24 = 18.4. We subtracted the nebular 
emission surrounding the object from the measured fluxes. 
The object HBHA 3703-01 was observed also during the 
2MASS survey (its catalogue ID is 20170810+3859295); its 
JHK S magnitudes are, respectively, J — 11.476 db 0.023, 
H = 10.458 ± 0.019, K s = 9.447 ± 0.017. The large colour 
index, J — K s — 2.03, tells that the star has an excess of 
IR emission most likely associated with a hot (>1000 K) 
circumstellar gas-dust envelope. 

In addition, observations of HBHA 370 3-01 were ob- 
tained by the AKARI infrared satellite (Ka wada et al.l 
12007m . According to these data, the fluxes from HBHA 3703- 
01 in the 65, 90, 140 and 160-/xm bands (in units of 
10" 16 erg cm-V 1 !" 1 ) are F 65 = 68.6, F 90 = 37.2, F140 = 
27.0 and Fi 60 = 18.3. Integration of SED in the 65-160- 
jLim range indicates that the total emission of dust in this 
range exceeds the luminosity of the central star by a factor 
of a few. Such a large flux may be due to the low angular 
resolution of AKARI (1-2 arcmin) JDoi Y. et alJl2012h . In 
addition to the circumstellar envelope, its field of view covers 
an appreciable part of the dust cloud visible on the 24-/xm 
image and heated by the radiation field of the Cyg OBI 
stars. Thus, the AKARI fluxes refer not only to the emis- 
sion of the circumstellar envelope of HBHA 3703-01 but also 
to the nebula GM 2-39 and to some part of its environment. 
In what follows we do not consider these fluxes. 

Figure E shows the SED of HBHA 3703-01 in the 0.44- 
24-/im spectral range plotted together with our photomet- 
ric data and 2MASS (J, if, K s ), as well as fluxes from 
the Spitzer images (3.5-24 ^ m). The figure als o shows the 
8-21-/mi MSX mission data (JEgan et alj|2003h . All fluxes 
were corrected for interstellar extinction that corresponds 
to E(B — V) = .87 mag using the extinction law from 
iRieke &; Lebofskvl (|1985T ). This value is less than the one 
used above in the estimation of the spectral type of the 
star, since some part of the extinction toward the central 
source may be due to the circumstellar envelope. The SED 
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Figure 8. The dereddened SED of HBHA 3703-01 in the 0.44- 
24 jum range (squares), MSX data (triangle) and calculated en- 
ergy distribution (solid curve). 



of HBHA 3703-01 taking into account the circumstellar dust 
envelope was calculated using th e CSDUST3 software, which 
was described in detail bv lEgan et al.l (|1988r ). We have done 
the following assumptions: 

1. The dust envelope was considered to be spherically 
symmetric with inner and outer radii R lu and Rout, respec- 
tively. The dust density distribution in the envelope was set 
as n(r) oc r~ a . 

2. The spectral energy distribution of the B(5-6)V cen- 
tral star was assumed to be blackbody with temperature 
T e ff = 15,000 K and luminosity L* = HOL©. Precisely this 
luminosity is obtained with a distance of 1.5 kpc and with 
the observed bolometric flux. 

3. The dust envelope consists of amorphous silicate and 
graphite particles with radius a S ii, a C ar and number densities 
risii and n ca r, respectively. The factors of efficiency of ab- 
sorption and scattering were calculated in accordance with 
the theory of Mie (1908) for spherical particles. The refrac- 
tive index of th e mate rial of the particles was taken from 
iDorschner et al.l d.1995). In the calculations we determined 
the optical depth of the envelope at a wavelength of 0.55 /im 
(tv), Rin, Rout, o S ii, a car , n S ii, n c&r and the exponent a in the 
dust density distribution law. Figure [8] presents the calcu- 
lated SED of the object HBHA 3703-01, which was obtained 
under the above assumptions. 

The best agreement between the computed and ob- 
served SED is achieved for the following parameters: ry — 
1.1, Ri n — 480 R*, a s n = 0.1 /mi, a car = 0.25 /mi, 
n si i/n ca r = 32, n{r) oc 1/r 1 " 5 , T in = 1350 K. 

Thus, in the circumstellar dust envelope of HBHA 3703- 
01 silicate particles prevail. The silicate dust temperature at 
the inner boundary of the envelope is Xi n =1350 K. Expo- 
nent a — 1.5 in the dust density distribution law is a typical 
value in accreti ng envelopes of Herbig Ae/Be stars (see, e.g., 
IShu et al.lll987T ) . The total mass of the dust envelope radiat- 
ing in the considered spectral range is Mdust ^3x 1O _6 M0. 
For the ratio M gas /Mdust ~ 100 this yields the total mass of 
the circumstellar envelope M env ~3x 1O _4 M0. Note that 
the estimated mass of the envelope depends on the outer ra- 
dius of the envelope; finding the envelope mass requires far- 



infrared observations with a high angular resolution. There- 
fore, we estimate only the mass of the material radiating in 
the wavelength range A < 25 /mi. 

Thus, the object HBHA 3703-01 located at the center 
of the diffuse nebula GM 2-39, in the region of the southern 
head of the IR globule, is a B(5-6)Ve star with a dust en- 
velope, whose temperature at the inner boundary is about 
1400 K. The star illuminates the nebula GM 2-39. The neb- 
ula GM 2-39 is distinguished from the ambient background 
of the molecular cloud just because of the presence of the 
above-mentioned star. 

Thus, the star HBHA 3703- 01 satisfies c ompletely the 
criteria for Herbig Ae/Be stars ()Herbig| [l960): (1) spectral 
type of A or earlier with emission lines, (2) location in an 
'obscured region' and (3) the star would illuminate a bright 
reflection nebulosity in its vicinity. 

In the area studied, no strong centimetre- wavelength 
radio continuum source is present that could be associ- 
ated with a compact ionized region. Only one compact ra- 
dio source has been detected there, MITG J201703+3900 
([Griffith et al.l Il99lh with flux density S„ = 53 mJy at 
v — 5 GHz. Most probably, it is a nonthermal extragalac- 
tic source, its spectral index a (S u oc v a ) in the 1.4-5 GHz 
interval is —2.32. 



6 CONCLUSION 

Our study of the interstellar medium toward a family of 
cometary globules detected in the archival data of the 
Spitzer Space Observatory in the region of a supershell in 
Cygnus allow us to draw the following conclusions. 

The analysis of kinematics of the interstellar medium 
in the neighbourhood of the family of CGs and of the stel- 
lar population in the regions suggests a physical association 
between CGs and the Cyg OBI association at a distance of 
1.5 kpc. (Ionized hydrogen and CO molecular clouds near 
the family of CGs have velocity Vlsr — 3-7 km s _1 ; the 
median radial velocity of the stars of Cyg OBI is Vlsr — 
4kms _1 .) 

The considered family of globules can be a remnant of a 
pillar directed toward the centre of the Cyg OBI association. 
The pillar is decaying in the process of star formation under 
the action of ionizing radiation and wind of the stars of 
Cyg OBI. 

The family of CGs is within the zone covered by the 
wind and ionizing radiation of the Of stars HD 228841 and 
HD 193595 belonging to Cyg OBI. The chain of CGs extends 
radially with respect to these stars. The morphology of CGs 
in the IR bands is as a whole consistent with the optical 
images: dust head and tail bordered with bright IR and/or 
optical rims. The optical emission encircles the bright IR 
globules. 

The ambient ionized gas is characterized by high- 
velocity features in the Ha line (up to 60-80 km s _1 with 
respect to the systematic velocity); this is a consequence of 
the action of the wind from these stars. The gas swept out by 
the wind circumnows the family of cometary globules, since 
high velocities are observed along the entire chain of CGs. 
We have found no obvious kinematic signs of braking of the 
high-velocity gas at the interaction with the 'Three-headed 
globule'. Indirect indications to the interaction of the wind 
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from the star and the family of CGs can be the presence of 
a CO cavern surrounding the eastern globules and fine IR 
filaments curled away from the stars. 

We have observed a compact stellar object in the south- 
ern head of the brightest globule - the star HBHA 3703-01 
identified with the infrared source IRAS 20153+3850. Our 
multicolour photometric observations allow us to estimate 
the spectral type of HBHA 3703-01 as B(5-6)V and to mea- 
sure its colour excess E(B — V) = 1.18 mag. Low-resolution 
spectral observations of HBHA 3703-01 have revealed the 
broad Ha emission with FWHM = 3.7 ± 0.5 A, which cor- 
responds to a high velocity, Av — 169 ± 25 km s _1 . The 
absolute visual magnitude of the star obtained under the 
assumption that the star belongs to the Cyg OBI starburst 
complex locates it on the H-R diagram near ZAMS in the 
region of Herbig Ae/Be objects. 

HBHA 3703-01 is surrounded by the reflection nebula 
GM 2-39; its diameter is about 30". The nebula is dis- 
tinguished against the background of the overall molecular 
cloud just because of the presence of the star, which illumi- 
nates a part of the interstellar medium nearest to the star 
and results in the formation of an optical reflection nebula. 

Using archival data of observations of HBHA 3703-01 
obtained by the Spitzer Space Telescope, we have estimated 
the fluxes in the 3.6, 4.5, 5.6, 8 and 24 /im bands. (On all IR 
images HBHA 3703-01 looks as a starlike object surrounded 
with a nebula; the emission of this nebula is subtracted from 
the measured fluxes.) We have also used our photometric 
data for HBHA 3703-01 together with the 2MASS magni- 
tudes (J, H, K 8 ). According to these data, we have derived 
the SED of HBHA 3703-01 in the 0.44-24 /xm range. All ob- 
servations were corrected for the interstellar extinction with 
E(B -V) = 0.87 mag. 

The obtained SED can be rather well described within 
the framework of the model of emission of the central source 
(a blackbody with temperature T = 15, 000 K) surrounded 
by a dust envelope with parameters listed above. 

The derived dust density distribution in the envelope 
n(r) oc 1/r 1 ' 5 is typical of the envelopes of Herbig Ae/Be 
stars. 



ACKNOWLEDGMENTS 

This work was based on the data of IFP observations in 
the Ha line, spectral observations on the 125-cm reflec- 
tor and photometric observations on the 60-cm reflector 
of the Crimean Laboratory of the Sternberg Astronomi- 
cal Institute, Moscow State University. We used archival 
data of the Spitzer, AKARI and 2MASS observations, ra- 
dio data in the CO line and the database of the Centre des 
donnees astronomiques, Strasbourg, France. This research 
was supported by the Russian Foundation for Basic Re- 
search (project codes 10-02-00091-a, 12-02-31356). OVE is 
also grateful for the financial support of the 'Dynasty' Foun- 
dation. We thank S.V. Antipin who made some of the pho- 
tometric observations. We thank the anonymous referee for 
the constructive comments that improved the quality of the 
paper. 



REFERENCES 

Bessell M.S., 1979, PASP, 91, 589 

Blaha C, Humphreys R. 1989, AJ, 98, 1598 

Brand P.W.J.L., 1981, MNRAS, 197, 217 

Brand P.W., Zealey W.J., 1977, A&A, 38, 363 

Brand P.W.J.L., Hawarden T.G., Longmore A.J., Williams 

P.M., Caldwell J.A.R., 1983, MNRAS, 203, 215 
Doi Y. et al., PKAS, 27, 111 
Dorschner J., Begemann B., Henning Th., Jager C, 

Mutschke H., 1995, A&A, 300, 503 
Droege T.F., Richmond M.W., Sallman M., 2007, PASP, 

119, 1083 
Egan M.P., Leung CM., Spagna G.F., 1988, Comp. Phys. 

Comm, 48, 271 
Egan M. P. et al., 2003, Air Force Research Laboratory 

Technical Report AFRL-VS-TR-2003-1589, VizieR online 

data catalogue, V/114 
Esipov V.F., Lozinskaya, T.A., Mel'nikov V.V., Pravdikova 

V.V., Sitnik T.G., Nichols-Bohlin J., 1996, Astron. Lett., 

22, 509 
Forbes D., English D., De Robertis M. M., Dawson P.C., 

1992, AJ, 103, 916 
Garmany CD., Stencel CD., 1992, A&AS, 94, 214 
Glushneva I.N., Doroshenko V.T., Fetisova T.S., Khruzina 

T.S., Kolotilov E.A., Mossakovskaya L.V., Ovchinnikov 

S.L., Voloshina LB., 1998, VizieR Online Data Catalog 

111/20^3 
Gottschalk M., Kothes R., Matthews H.E., Landecker T.L., 

Dent W.R.F., 2012, A&A, 541, A79 
Griffith M., Heflin M., Conner S., Burke B., Langston C, 

1991, ApJS, 75, 801 
Gyulbudaghian A.L., Maghakian T.Y., 1977, Doklady 

Akad. Nauk Armenian SSR, 64, 104 
Hawarden T.G., Brand P.W.J.L., 1976, MNRAS, 175, 19 
Herbig G.H, 1960, ApJS, 4, 337 
Hora J.L. et al., 2008, New Light on Young Stars: Spitzer 's 

View of Circumstellar Disks, Proc. Fifth Spitzer Conf. 

CalTech, Pasadena, CA0 
Kawada M. et al., 2007, PASJ, 59, S389 
Kniazev A.Y., Pustilnik S.A., Zucker D.B., 2008, MNRAS, 

384, 1045 
Kohoutek L., Wehmeyer R., 1999, A&AS, 134, 255 
Landecker T.L., 1983, JRAS Canada, 77, 264 
Landolt A., 1992, AJ, 97, 337 
Leung H.O., Thaddeus P., 1992, ApJS, 81, 267 
Lozinskaya T.A., Sitnik T.G., 1988, SvA Lett., 14, 100 
Lozinskaya T.A., Repin S.V., 1990, SvA, 34, 580 
Lozinskaya T.A., Pravdikova V.V., Sitnik T.G., Esipov 

V.F., Mel'nikov V.V., 1998, Astron. Reports, 42, 453 
Mie C, 1908, Ann. Physik, 25, 377 
Pittard J.M., 2011, MNRAS, 411, L41 
Reipurth B., 1983, A&A, 117, 183 
Richmond M.W. et al., 1994, A J, 107, 1022 
Rieke G.H., Lebofsky M.J., 1985, ApJ, 288, 618 
Sahu M., Pottasch S.R., Sahu K.C, Wesselius P.R., Desai 

J.N, 1988, A&A, 195, 269 
Sandquist A, 1976, MNRAS, 177, 69P 



5 http : //cdsarc . u-strasbg . f r/viz-bin/Cat?cat=III°/ 2F208& 

6 http : //www . ipac . caltech . edu/spitzer2008/ 
talks/ JosephHora_Spitzer2008Conf erence . pdf 



12 V. P. Arkhipova et al. 

Sergeev S.G., Heisberger F., 1993, A Users Manual for SPE. 
Wien 

Shu F.H., Adams F.C., Lizano S., 1987, Annual Rev. As- 
tron. Astrophys., 25, 23 

Sitnik T.G., 2003, Astron. Lett., 29, 311 

Sitnik T.G., 2011, Astron. Reports, 55, 616 

Sitnik T.G., Lozinskaya T.A., 2009, Astron. Lett., 35, 121 

Sitnik T.G., Mel'nik A.M., 1996, Astron. Lett., 22, 122 

Smith L.J., Norris R.F. Crowther P.A., 2002, MNRAS, 337, 
1309 

Smith N., Bally J., Walborn N.R., 2010a, MNRAS, 405, 
1153 

Smith N. et al., 2010b, MNRAS, 406, 952 

Stetson P.B., 2000, PASP, 112, 925 

Straizys V., 1977, Multicolor Stellar Photometry. Photo- 
metric Systems and Methods. Mokslas Publishers, Vilnius 

Straizys V., 1982, Metal-deficient Stars. Mokslas Publish- 
ers, Vilnius 

Turner D.G., 1989, 1989, AJ, 98, 2300 



This paper has been typeset from a TgX/ WF^K. file prepared 
by the author. 



